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1.0  -  INTRODUCTION 

This  report  provides  a  methodology  for  determining  analyticaUy  the  Receiver  Operating 
Characteristic  (ROC)  curves  of  the  Mid-Frequency  Active  Classification  Processor  (MFACP) 
Continuous  Wave  (CW)  normalizer.  The  methodology  is  based  on  the  determination  of  the 
statistics  of  the  normalizer  estimate  under  stationary  and  non-stationary  noise  conditions. 

The  model  assumes  the  averaging  time  of  the  normalizer  is  long  enough  so  that  the 
probability  density  function  of  normalizer  estimate  p  converges  to  a  Gaussian  form  as  a 
consequence  of  the  Central  Limit  Theorem.  The  dynamic  noise  background  model  described  in 
section  4.0  has  a  sinusoidal  variation  in  the  Rayleigh  parameter  of  the  envelope  deteaed 
narrowband  Gaussian  noise.  This  model  of  the  shape  of  the  background  is  quite  general  and  in 
reality  the  shape  of  the  background  variation  could  be  arbitrarily  specified  as  required  as  long  as 
the  requirements  given  in  4.0  are  met.  Section  2.0  describes  the  false  alarm  and  detection 
probabilities  for  each  test  cell  under  investigation  which  leads  to  the  generation  of  the  ROC 
curves.  Section  3.0  develops  the  model  for  determining  the  mean  and  variance  of  the  normalizer 
estimate  and  section  4.0  describes  the  background  noise  model. 
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2.0  -  PERFORMANCE  EVALUATION  METHODOLOGY 

The  perfoimance  of  each  nonnalization  algorithm  is  determined  by  evaluating  the 
Receiver  Operating  Characteristic  (ROC)  curves.  These  curves  quantify  the  performance  of  the 
algorithm  as  a  function  of  Signal-to-Noise  Ratio  (SNR)  and  probability  of  detection,  P(D),  at 
different  probabilities  of  false  alarm,  P(F). 

A  fimdamental  requirement  of  a  normalization  algorithm  is  a  Constant  False  Alarm  Rate 
output,  CFAR.  The  MFACP  nuiintains  a  constant  detection  threshold  at  every  range  point  in  the 
detection  process.  The  CFAR  requirement  is  achieved  when  the  noise  background  is  stationary, 
i.e.,  the  statistics  do  not  change  with  time.  In  non-stationary  noise  backgrounds  with  Exed 
detection  thresholds  the  single  bin  P(F)  and  P(D)  is  dynamic  and  therefore  the  CFAR 
requirement  is  not  met. 

The  test  cell  x  is  divided  by  the  background  mean  level  estimate  p  generating  the 
normalized  output.  Let  X  denote  the  deteaion  threshold.  Ho  the  hypothesis  that  the  test  cell  x 
contains  noise  ordy  and  Hi  the  hypothesis  that  test  cell  x  contains  a  target  echo  (signal  plus 
noise).  (For  notadonal  convenience  the  time  dependence  for  each  successive  test  cell  has  been 
dropped).  The  following  test  is  carried  out  to  determine  a  detection: 


< 


P(F)  is  the  declaration  that  a  true  detection  was  made  under  hypothesis  Hq,  i.e.,  a  detection  is 
erroneously  declared  since  x  contained  noise  otily.  P(D)  is  the  declaration  that  a  true  detection 
was  made  under  hypothesis  Hi,  i.e.,  a  detection  is  correctly  declared  since  x  contains  signal  plus 
noise.  Under  the  noise  only  hypothesis  Hq,  the  probability  density  function  of  the  envelope  of  a 
narrowband  Gaussian  process  has  a  Rayleigh  form  given  in  [7]  as 


f(x  I  Hq)  =  ^  exp  u(x) 
o  2a 


u(x)  is  the  unit  step  function.  For  a  narrowband  signal  in  narrowband  noise,  i.e.  signal  plus 
noise  hypothesis  Hi,  the  probability  density  function  of  the  test  cell  x  is  Rician  and  given  by 


f(x  I  Hi) 


(2-3) 
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where  1©  (•)  is  the  modified  Bessel  function  of  order  zero.  The  SNR  in  dB  for  the  test  cell  x 
containing  a  target  echo  is 


SNR  =  10  logjo 


r  a2\ 


(2-4) 


It  is  assumed  that  the  normalizer  window  size  (averaging  time)  is  large  enough  so  that  the 
convergence  property  of  the  Central  Limit  Theorem  holds  and  that  a  Gaussian  model  for  the 
normalizer  estimate  p  is  valid.  The  probability  density  funaion  of  the  normalizer  estimate  p  is 
[4-6] 


1/V^q^ 

FG(mji/Oji) 


/ 

exp 


U(p) 


(2-5) 


where  Fg(»)  is  the  cumulative  distribution  of  the  standard  normal  zero  mean,  unit  variance  form 
and 


m^  =  E[p] 
aj  =  var(p). 


P(F)  and  P(D)  are  also  given  in  [4-6]  as 


P(F)  = 


^  XI  Ho 


'  J  1 

Fg 

f  y 

^^H-X^aJ/a^  J 

“‘t  2  l+XV./O^J 

FG(mji/Oj^) 

(.2-6) 


(2-7) 


and 
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P(D)  =  ^'HiJ 

=  £2Q(A  /  c,  /  o)  dp.  (2-8) 


where  Q(qo,  qi)  is  the  Marqum-Q  function  for  arguments  qo  and  qi  defined  by 

Q(qo.  qi)  =  j^Vexp(-(V^+qo)/2)Io(qoV)dV 


(2-9) 


Wn  and  Zn  are  the  Nh  weights  and  zeros  respectively  associated  with  the  Hermite  polynomial 
expansion  and  are  well  tabulated  in  [9]  for  Nh  up  to  20.  The  polynomial  expansion  up  to  Nh  =  6 
provides  very  high  accuracy.  The  Rayleigh  parameter  a  given  in  equations  (2-7)  and  (2-8)  is  that 
of  the  test  cell  under  investigation. 

The  analytic  determination  of  the  time  history  of  the  noimalizer  output  mean  and 
variance  which  are  necessary  in  the  determination  of  the  RCK)  curves  is  given  in  section  3.0. 
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3.0  -  MFACP/CW  NORMAUZER  FUNCTIONAL  DESCRIPTION 

Figure  (3-1)  is  a  block  diagram  of  the  signal  processing  prior  to  normalization.  It  consists 
of  (a)  a  windowing  fimction  which  is  either  a  rectangle  or  Hanning,  (b)  a  Fast  Fourier  Transform 
(FFT)  with  75%  overitq)  and  (c)  an  envelope  deteaor. 

Figure  (3-2)  is  a  block  diagram  of  the  MFACP/CW  normalizer.  The  mean  level  at  each 
range  bin  for  a  given  doppler  chaimel  is  estimated  by  an  exponential  filter.  The  min  function 
(minimum  of  two  values)  is  used  for  large  data  outlier  rejection  so  that  the  mean  background 
level  estimate  is  not  biased.  The  given  clipping  threshold  Xc  of  2.56  provides  an  optimal  rate 
limitation  [2, 3]  of  0.76  dB  per  FFT  update. 

The  mean  background  level  estimate  in  each  doppler  channel  at  time  k,  denoted  by  z(k), 
is  the  sum  of  the  scaled  input  at  time  k,  x(k),  and  die  scaled  mean  background  level  output  at 
time  k-1,  z(k-l).  The  input  data  to  the  exponential  filter,  x(k),  is  the  minimum  of  the  envelope 
detected  output  at  time  k,  x(k),  and  the  mean  background  level  estimate  z(k-l).  The  exponential 
filter  recursively  updates  the  mean  background  level  estimate  at  each  time  sample  k  using  the 
clipped  data  at  the  output  of  the  minimum  function. 

To  further  reduce  the  variance  of  the  mean  background  level  estimate,  5  doppler  channels 
arc  averaged  before  generating  the  final  normalizer  estimate.  The  average  of  the  five  doppler 
channels,  z'(k),  is  symmetric  about  the  channel  of  interest;  i.e.  2  adjacent  channels  on  either  side 
of  the  channel  of  interest  are  used.  To  nullify  the  bias  in  the  mean  background  estimate  z'(k) 
when  the  test  cell  of  interest  contains  a  target,  a  range  gap  symmetric  about  and  centered  on  the 
test  cell  is  normally  used.  The  MFACP/CW  normalizer  is  one-sided  (asymmetric)  and 
implements  the  gap  by  normalizing  the  test  cell  with  the  mean  background  estimate  delayed  by  4 
samples,  z'(k-4). 

In  the  analysis  that  follows  the  mean  and  variance  of  the  output,  z'(k),  are  generated.  The 
dynamic  behavior  of  these  two  values  determine  the  ROC  curves  at  each  range  bin  in 
nonstationary  noise  environments. 
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Figure  (3-1)  -  Signal  processing  prior  to  normalization. 


3.0.1  -  MFACP/CW  CORRELATION  COEFnCIENTS 


For  overli^ped  FFT  processing  as  shown  in  figure  (3-1)  above  the  correlation  coefficients 
between  successive  updates  [1]  is 


X^!I,‘‘w(n)  W(n  +  (l-n)N) 


P(n)  = 


(3-1) 


O  is  the  fractional  overlap  and  W(n)  is  the  window  function.  Figure  (3-3)  shows  the  sequencing 
of  the  data  through  the  FFT  for  75%  overltq).  After  4  updates  the  new  input  sequence  to  the  FFT 
has  no  data  points  in  common  with  the  original  input  sequence. 

For  the  four  updates  die  fractional  overly  O  is  0.75, 0.5,  0.25  and  0  respectively  for  i  = 
1, 2, 3, 4.  The  correlation  coefficient  at  the  envelope  detector  output  is  given  in  [5, 6]  as 


r(i) 


P^i)/4 
4/7C  -  1 


n  + 


P^(i)y«”» 

4  “^n=0 


2"(n-H)! 


P^“(i)[ 


(3-2) 


where  i  =  1, 2,  3, 4  corresponds  to  the  FFT  update.  The  infinite  sum  in  equation  (3-2)  may  be 
truncated  to  a  small  finite  number,  Nn,  of  terms  producing  great  accuracy  (Nn  S  30).  Table  (3- 
1)  lists  the  correlation  coefficients  at  the  outputs  of  the  FFT  and  envelope  detector  for  both  the 
rectangular  and  Hanning  windows.  Figure  (3-4)  is  a  plot  of  the  correlation  coefficients 
comparing  both  window  types  at  the  FFT  and  envelope  detector  outputs. 
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Update,  i 

pp,  rectangle 

r^  rectanlge 

p„  Hanning 

r^  Hanning 

0 

1.0 

1.0 

1.0 

1.0 

1 

0.75 

0.5361 

0.6592 

0.4098 

2 

0.50 

0.2326 

0.1667 

0.0255 

3 

0.25 

0.0574 

0.0075 

0.0001 

4 

0 

0 

0 

0 

Table  (3-1)  -  Conelation  coeffideots  at  output  of  FFT  and  envelope  detector. 


The  noise  background  at  the  FFT  output  is  assumed  to  be  a  zero  mean  Gaussian  process,  N(0, 
a2).  At  the  output  of  the  envelope  detector  the  noise  is  a  Rayleigh  process  [7]  with  Rayleigh 
parameter  O"  equal  to  the  variance  of  the  Gaussian  noise.  The  density  function  for  a  Rayleigh 
random  variable  is 


/ 

p(x)=  exp 


-X 


2  ^ 


u(x) 


(3-3) 


where  u(x)  is  the  unit  step  function.  Figure  (3-5)  is  an  example  plot  of  equation  (3-3)  for  the 
Rayleigh  parameter  =  1.  The  mean  and  variance  of  a  Rayleigh  random  variable  x  are 

=  ^|n/2  a 
oj  =  (2  -  7c/2)<t^. 


original  ii^ut  segment 

=  0.75  update  1 

£2  =  0.50  update  2 

Q  =  0.25  update  3 

f2  =  0  update  4 

Hguie  (3-3)  -  Sequencing  of  data  segments  for  75%  overbq)  FFT  processing. 


correlation  coefficient 
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Figure  (3-4)  •  Plot  of  the  correlatioa  coefficients  of  FFT  update  at  the  output  of  the  oveilapped  FFT 
(equatitm  Al)  and  envelope  detector  (equation  3-2)  for  a  rectangular  and  Hanning  window  as  given 
in  table  (3-1). 


Figure  (3-S)  -  Example  plot  of  the  Rayleigh  density  function  of  the  samples  at  the  output  of  the 
envelope  detector  which  are  input  to  minimum  function  for  Rayleigh  parameter  s  i. 
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3.0.2  -  EXPONENTIAL  FILTER 

The  exponential  filter  [2, 3]  ouqiut  at  time  k  (see  Hgure  (3-2))  denoted  by  z(k)  is 
z(k)  =  ax(k)  +  pz(k-l)  (3^) 

where  a  =  — ^ — ,  B  =  ^  -  andN  =  15. 

N  +  1  ^  N  +  1 

Modifying  z(k)  to  incoiporate  M+1  previous  input  values  gives 

+  P^*^z(k-M-1).  (3-5) 

The  mean  of  the  exponential  filter  output  z(k)  denoted  by  ^^(k)  =  E[z(k)]  is 

p^(k)  =  E[z(k)]  =  aE[x(k)]  +  pE[z(k.l)] 

=  apx(k)  +  P)i2(k-1)  (3-6) 

=  +  P’^*^E(z(k-M.l)]. 

The  variance  at  the  exponential  filter  output  denoted  by  a^(k)  is  given  by 
a?(k)  =  E[z^k)]  -  E[z(k)f 

2  2  o2  2  « 

=  a^aJCk)  +  p^a‘(k-l)  +  2apcov(x(k),  z(k-l)). 

The  covariance  between  the  present  input  x(k)  and  the  previous  output  z(k-l),  denoted  by 
cov(x(k),  z(k-l))  in  equation  (3-7),  is  by  definition  [8] 

cov(x(k),z(k-l))  =  E[x(k)  z(k  - 1)]  -  E[x(k)]  E(z(k  - 1)]. 

Substituting  the  relationship  given  in  equation  (3-5)  for  the  previous  M+1  input  values  the 
covariance  expression  becomes 
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cov(x(k).  z(k-l))  =  (E[x(k)i(k-l-i)]  -  E[x(k)]E[x(k-l-i)]) 


(3-8) 


where  also  by  definition  [8] 

cov(x(k).  x(k-l-i))  =  r(i)Oj(k)Oi(k-l-i) 

was  used.  Substituting  equation  (3-8)  back  into  (3-7)  and  the  faa  that  the  correlation  coefficient 
r(i)  goes  to  zero  after  M  +  1  =  4  updates  the  etqxtnential  filter  output  variance  becomes 

a^(k)  =  a^<sl(k)  +  p2aj(k-l)  -h  2a^^“^^P*ajt(k)Ojj(k-i)r(i).  (3-9) 

Equations  (3-6)  and  (3-9)  are  general  expressioiu^  ^  ‘  'mg  the  recursive  relationship  for  the  mean 
and  variance  of  the  exponential  filter  output  in  either  stationary  or  non-stationary  noise 
backgrounds.  If  the  filter  is  operating  in  a  stationary  noise  field,  the  filter  goes  to  steady  state 
and  the  output  mean  and  variance  becomes 


^2(k)  =  lix(k)  =  \i^ 
a^ai 


o^(k) 


1-P' 


(l  +  2  = 


(3-10) 


To  reduce  the  variance  of  the  exponential  filter  ouqjut  5  doppler  channels  are  averaged.  The 
variance  is  thus  reduced  to,  by  assuming  the  5  doppler  chaiuiels  are  independent  and  identically 
distributed, 

alik)  =  a2(k)/5.  (3-11) 


The  input  test  cell  under  investigation  is  divided  by  the  mean  estimate  delayed  by  4  samples 
producing  the  normalized  ouq>ut 


y(k)  = 


x(k) 

z’(k-4)‘ 


(3-12) 
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3.0.3  -  MFACP/CW  MINIMUM  FUNCTION 


The  minimum  function  chooses  the  smaller  of  the  Rayleigh  output  from  the  envelope 
deteaor  at  time  k  and  the  scaled  output  of  the  exponential  filter  at  time  k-1 .  This  function 
essentially  shears  (clips)  the  data  to  eliminate  large  data  outliers  from  corrupting  the  mean  level 
estimate  and  provides  a  rate  limitation  of  0.76  dB  per  FFT  update  [2, 3].  The  density  function  of 
a  clipped  Rayleigh  random  variable  is 


Pj(x)  =  -y  c 

c 


(  (  _y2  \ 

^  TT  1“^*)  *  +  exp 

\^2a  j  yia  j 


(3-13) 


where  5(x  -  )  is  the  unit  impulse  funaion  located  at  x  =  amplitude  of  the  delta 

function  is  the  probability  of  the  Rayleigh  random  variable  exceeding  the  clipping  threshold. 
The  mean  and  variance  of  the  clipped  Rayleigh  random  variable  are  given  by 


Vi.  =  E[xl 


■© 


-I  oerf 


E[x^]  =  2a^  1  -  ex] 


P  — T 


(3-14) 


<4  =  E[x2]  -  E[x]2. 


erf(<t>)  is  the  standard  error  function  defined  in  [9]  as 


erf(<J>)  =  J^exp(-Y^)dY. 


(3-15) 


A  plot  of  the  clipped  Rayleigh  density  function  specified  by  equation  (3-13)  is  shown  in  figure 
(3-5)  for  =  1  and  =  2.  The  height  of  the  delta  fimction  is  equal  to  the  probability  of  the 
Rayleigh  random  variable  x  >  Xc. 

3.0.4  -  MFACP/CW  EXPONENTIAL  FILTER  OUTPUT  NOISE  MODEL 


The  exponential  filter  output  mean  and  variance  is  determined  by  combining  the  results 
from  3.0.2  and  3.0.3.  To  determine  the  mean  and  variance  of  the  exponential 
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Hgure  (3-6)  •  Example  plot  of  the  density  funcdon  of  the  eloped  Rayleigh  random  variables  at  the 
on^ut  of  the  minimum  function  for  Rayleigh  parameter  =  1  and  Xc  =  2.  The  amplitude  of  the 
delta  function  is  equal  to  the  probability  of  the  Rayleigh  random  variable  x  ^  Xc . 


output  the  algorithm  as  shown  in  figure  (3-2)  is  combined  with  equations  (3-6, 3-9)  and  (3- 
14).  This  model  results  in  the  following  expressions  for  the  mean  and  variance  of  the 
exponential  filter  output  at  each  time  k: 


Pj(k)  =  amin(p.3^(k),  2.56n2(k-l))  +  pp^fk-l) 
=  m,, 


(3-16) 


o|(k)  =  2<rJ(k)  |l  - 


.2\ 


-(fp.(k)(erf 


:  -D) 

(k) 

2.56p^(k-l)'|f 

,  V2o,(k)  I 


yj 


(3-17) 


a?(k)  = 


=  o 


(aa^(k)f  +  (pa^(k-l))2  + 
2 


(?-i8) 
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3.0J  -  EFFECTIVE  AVERAGING  TIME 


The  effective  noimalizer  averaging  time  is  determined  as  the  number  of  independent 
samples  that  go  into  the  computation  of  the  mean  estimate  for  random  error  reduction.  For 
stationary  noise  with  c  =  1  the  mean  and  variance  of  the  clipped  Rayleigh  random  variable  using 
equation  (3*14)  evaluates  to 


\L.  =  E[x]  = 


2.56Vir/2 


a?  =  2 


l-exp|^ 


-2.56^(7c/2) 


V2 

X 


1.25163 


1.25163 

^  =  0.42178. 


(3-19) 


Substituting  now  into  equation  (3-10)  gives  the  variance  of  the  estimate  before  doppler  averaging 
as 


0.42178 

31 


(1.81336) 


0.42178 

17 


(3-20) 


The  effective  averaging  time  of  the  recursive  stracture  of  the  normalizer  is  deteimined  by 
comparing  to  an  equivalent  window  structure  as  would  occur  in  a  split  window  scheme.  Looking 
at  equation  (3-11)  the  equivalent  length  of  the  split  window  block  averager  is  given  by 


Nno  = 


(I-P^)Ndoppler  ^  (31X5) 
a^(l  +  2]£*^J’^P‘r(i))  1-81336 


(3-21) 


where  Ndoppler  is  the  number  of  Doppler  channels  averaged.  Therfeore  by  averaging  the  5 
doppler  channels  the  effective  averaging  time  is  increased  from  Nno  =  17  to  85  samples. 


3-10 


TM  921205 


4.0  -  BACKGROUND  NOISE  MODEL 


The  noise  nonstationaiity  considered  here  [5]  allows  for  an  arbitrary  positive  sample-to- 
sample  variation  in  the  Rayleigh  parameter  Ok-  The  mean  of  the  single  noise  cell  varies  from 
sample-to-sample  as 

=  E[Xk]  =  Vn/2  0k  (^1) 


This  type  of  nonstationarity  corresponds  to  a  sample-to-sample  change  in  the  total  noise  power 
without  necessarily  a  variation  in  the  power  spectrum  with  frequency.  The  latter  nonstationarity, 
involving  a  spectral  shape  change  without  a  total  power  change,  does  not  afreet  the  first  order 
statistic,  E[xk],  but  does  affect  the  sample-to-sample  covariance 

Vij  =  cov(Xi,Xj)  =  rjj  (4-2) 

where  is  the  variance  of  sample  Xk  where  k  =  i,  j 

aj.  =  (2  -  7c/2)of  (4-3) 


by  changing  the  correlation  coefficient,  tij. 


"ij  = 


E[XiXj]-E[xJE[Xj] 


(4^) 


this  total  power  nonstationarity  also  affects  Vy  through  and  .  Thus,  noise  nonstationarity 
is  modeled  by  a  total  power  nonstationarity  afrecting  both  m^.  and  v^j,  and  a  spectral  sh^ 

nonstationarity  with  frequency  affecting  only  Vy.  In  stationary  noise  backgrounds,  the 
correlation  coefficients  ty  does  not  vary  with  time.  In  nonstationary  noise,  the  Rayleigh 
parameter  Oi  change  from  sample-to-sample  and  the  correlation  coefficients  may  change  as  well. 

In  shallow  water  enviromnents  it  has  been  observed  in  real  data  that  the  [6]  background 
variation  is  oscillatory  in  nature  and  can  vary  by  as  much  as  20  dB  in  the  space  of  a  few  seconds. 
In  convergence  zone  regions,  the  background  may  vary  by  as  much  as  15  to  20  db  per  second 
and  has  a  pedestal  shape  to  the  variation.  The  effects  on  normalizer  performance  is  captured  by 
considering  a  generalized  sinusoidal  background  variation  model  where  both  the  amplitude  and 
period  of  the  variation  may  be  arbitrarily  specified. 
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Let  the  variation  in  the  Rayleigh  parameter  as  a  function  of  tiine  be  defined  by 

o(t)  =  Co  -  Omin  sin(2jct/T)  (4-5) 

where  d  is  the  total  power  change  in  dB  and  T  is  the  period  of  the  sinusoid.  Also  let  d  be  deHned 
by 

d  =  201ogio(o„„/o^)  (4-6) 


or 


=  o, 


max 


/  ^mir 


(4-7) 


so  that  at  the  minimum  and  maximum  points  in  the  sinusoidal  variation 

Ojoax  -  Oo-<yisin(37c/2)  =  ao+Ci 
<ymin  =  <Jo-<JlS>n(Jt/2)  =  Oo-Ci. 

Equating  equations  (4-7)  and  (4-8)  gives 


10<i/2O 

Co+Ci 

1 

Q 

Cl  - 

10^/20  ^.1^ 

and  finally  the  variation  of  the  Rayleigh  parameter  for  the  sinusoid  can  be  written  as 
c(t)  =  Co 


r 

flO'*'*’-!') 

1 

fH 

_ 1 

10<l/2O^l^ 

H-rJ 

(4-8) 


(4-9) 


(4-10) 
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